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The principal axes of the recently discovered anisotropic phases of 2D electron systems at high 
Landau level occupancy are consistently oriented relative to the crystal axes of the host semicon- 
ductor. The nature of the native rotational symmetry breaking field responsible for this preferential 
orientation remains unknown. Here we report on experiments designed to investigate the origin and 
magnitude of this symmetry breaking field. Our results suggest that neither micron-scale surface 
roughness features nor the precise symmetry of the quantum well potential confining the 2D system 
are important factors. By combining tilted field transport measurements with detailed self-consistent 
calculations we estimate that the native anisotropy energy, whatever its origin, is typically ~ 1 mK 
per electron. 



A fundamental goal of contemporary condensed mat- 
ter physics is to understand the ground state of the two- 
dimensional electron system (2DES) and to explore its 
experimental signatures. In the presence of a large per- 
pendicular magnetic field B, the kinetic energy of 2D 
electrons becomes quantized into discrete massively de- 
generate Landau levels (LLs) . At high B only the lowest 
(N = 0) LL is occupied, and the 2DES exhibits its most 
spectacular phenomenon: the fractional quantized Hall 
effect (FQHE) Q. At lower magnetic fields the higher 
LLs become occupied, and in the third or higher (N > 2) 
LL no FQHE states have been observed. Recent exper- 
iments have nevertheless uncovered extraordinary trans- 
port signatures unique to the high LL regime, pointing to 
the existence of a new class of many-body states distinct 
from the incompressible quantum fluids responsible for 
the FQHE |^||. The most striking of these signatures 
is the rapid development, at very low temperatures, of 
strong anisotropies in the longitudinal resistance of the 
2D electron system. 

Observed only in very high mobility samples and at 
temperatures below about 100 mK, the anisotropies in 
the longitudinal resistance are strongest near half-filling 
of the N — 2 and several higher Landau levels. This 
corresponds to Landau level filling fractions v, defined 
as the ratio of the electron density N s to the degener- 
acy eB/h of a single spin-resolved LL, of v = 9/2, 11/2, 
13/2, etc. For 2D electron systems in GaAs/AlGaAs het- 
erostructures grown on [001]-oriented GaAs substrates, 
the anisotropies are consistently disposed so that the 
"hard" transport direction is parallel to the [110] crys- 
tallographic direction while the "easy" direction is par- 
allel to [110]. Although a persuasive picture of how 
anisotropic electronic ground states develop at high Lan- 
dau level occupancy now exists, there is still no under- 
standing of why they are consistently oriented relative to 
the host crystal axes. The origin and magnitude of the 



necessary native symmetry breaking field is the focus of 
this paper. 

The resistance anisotropy near half-filling of high LLs 
has been widely interpreted as evidence for charge- 
density- wave (CDW) ground states. Hartree-Fock (HF) 
calculations |||| have suggested that 2D electrons in 
half-filled high LLs form a unidirectional CDW, or 
"stripe" phase. At v = 9/2, for example, the system 
is expected to phase separate into alternating regions of 
v = 4 and v — 5 with a period of about three times 
the cyclotron radius R c = Tikp/eB, or about 100 nm in 
typical samples. If these stripes are somehow preferen- 
tially oriented, anisotropy in the longitudinal resistance 
of the sample would likely result. Transport currents 
flowing perpendicular to such stripes would presumably 
encounter greater resistance than currents flowing paral- 
lel to them. 

More recent theoretical work has generally supported 
the early HF predictions, albeit with some significant 
modifications. Numerical exact diagonalizations of small 
systems of electrons have uncovered sharp features in the 
susceptibility and structure factor of the 2DES in high 
LLs consistent with stripe formation at half filling [|6|. 
Refinements to the HF approach suggest that the stripes 
are themselves unstable against modulations along their 
length [0J|]. The resulting "smectic crystal" would pre- 
sumably be pinned and therefore insulating at very low 
temperatures. Quantum and/or thermal fluctuations are 
expected to melt the system and render it analogous to 
a nematic liquid crystal ||. A nematic has translational 
but not rotational symmetry and might therefore exhibit 
anisotropic transport if a preferred direction exists. At 
still higher temperatures a nematic to isotropic transition 
is predicted to restore rotational symmetry ]9],|To|. 

The consistency of the orientation of the transport 
anisotropy in high LLs has been established through 
experiments on a large number of GaAs/AlGaAs het- 
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erostructure samples. Grown by molecular beam epi- 
taxy (MBE), these modulation-doped 2D electron sys- 
tems have low temperature mobilities typically around 
10 7 cm 2 /V s and densities ranging from 1.5 — 2.7 x 
10 11 cm~ 2 . The samples are usually illuminated briefly 
at low temperature with a red light emitting diode to en- 
hance their transport characteristics but, as explained 
below, the existence and orientation of the high Lan- 
dau level resistance anisotropy does not depend upon 
this. Although most measurements have been made 
using square samples whose sides are parallel to the 
[110] and [110] directions, specific additional experi- 
ments have shown that the orientation of the transport 
anisotropy is insensitive to the geometry of the conduct- 
ing region, at least in mm-sized samples jll] . 

In this paper we report new results which narrow down 
the list of possible native symmetry breaking effects that 
might orient the high LL anisotropy. Specifically, we 
show that neither the micron-scale morphology of the 
sample surface, as revealed by atomic force microscopy 
(AFM), nor the shape of the heterojunction potential 
confining the electron gas to two dimensions appear to 
be significant factors. In addition, by combining experi- 
mental results obtained using tilted magnetic fields with 
detailed numerical calculations, we estimate the strength 
of the still unknown symmetry breaking field. 

Data from three samples (A, B, and C) will be pre- 
sented here. Sample A consists of a single 30 nm 
GaAs quantum well embedded in thick layers of Al- 
GaAs. Doping is provided by Si layers positioned sym- 
metrically above and below the quantum well. Before 
illumination the 2DES in this sample has a density of 
2.74 x 10 11 cm -2 and a low temperature mobility of 
1.8 x 10 7 cm 2 /V s. Brief illumination at 1.6 K changes 
these values to 2.54x 10 11 cm -2 and 2.1 xlO 7 cm 2 /V s, re- 
spectively. Samples B and C are conventional single het- 
erointerfaces. After illumination these samples respec- 
tively have densities of 2.67 and 1.48 x 10 11 cm" 2 and 
mobilities of 0.9 and 1.1 x 10 7 cm 2 /V s. Each sample 
was cleaved from its parent MBE wafer into a square, 3 
or 5 mm on a side, and eight diffused In contacts were po- 
sitioned at the corners and midpoints of the sides. Trans- 
port measurements were made with an excitation current 
of / < 20 nA at 13 Hz driven between midpoint contacts 
on opposite sides of the sample. Longitudinal resistance 
data were obtained by measuring the resulting voltage 
drop between corner contacts. Each of these samples 
clearly exhibits the characteristic low temperature resis- 
tance anisotropics near half-filling of the N > 2 LLs first 
reported by Lilly, et al. M. 

The inevitable imperfections which arise during crystal 
growth are natural candidates for an extrinsic mechanism 
capable of orienting the anisotropic electronic phases in 
high LLs. Since the kinetics of MBE growth is known to 
be anisotropic, it is easy to imagine that these imperfec- 
tions might be preferentially oriented. This would break 



the rotational symmetry of the local environment of the 
2DES and possibly influence its transport properties. As 
it is obviously difficult to directly observe the local envi- 
ronment of the 2DES, which is buried inside the sample, 
we have resorted instead to an examination of the ex- 
posed top surface of the heterostructure via AFM Jl2| . 
Although this surface is typically 150 nm above the plane 
of the 2DES, the predicted wavelength of the CDW in 
high LLs is of comparable magnitude. 

Nine high mobility GaAs heterostructure wafers (in- 
cluding those of samples A, B, and C) were examined 
in search of anisotropic surface morphologies. Each con- 
tained a 2DES with mobility exceeding 0.8 x 10 7 cm 2 /V s 
and density in the range of 1.5 — 2.7 x 10 11 cm -2 . Unam- 
biguous determinations of the [110] and [110] directions 
were made using one or more techniques [|l3|. Samples 
from all nine wafers showed strong resistance anisotropics 
at v = 9/2, 11/2, etc. with the high resistance direction 
along [1T0]. 

Figure 1 shows the results of transport and AFM stud- 
ies of samples A and B. The top two panels show longitu- 
dinal resistance data, taken at T = 50 mK, from the fill- 
ing factor range 3 < v < 5. Near v — 9/2 the strong resis- 
tance anisotropy characteristic of half-filled high Landau 
levels is clearly evident. As the data show, the "hard" 
transport direction is along [110]. In contrast, the resis- 
tance near v = 7/2 is nearly isotropic. This filling factor, 
which lies in the N = 1 LL, supports a weak fractional 
quantized Hall state which gives way to an anisotropic 
state only when a large magnetic field component paral- 
lel to the 2D plane is applied [ll]j2"2] |. It is worth noting 
that the data from sample A shown here was obtained 
without illumination, verifying that the anisotropy effect 
does not depend upon this. Data from sample A after 
illumination will be discussed below. 

The lower panels of Fig. 1 display 16 x 16 /im AFM 
images of the [001] -oriented surfaces of the wafers from 
which samples A and B were cleaved, and the [110] and 
[110] directions arc indicated by arrows. The images were 
obtained at room temperature using the "tapping mode" 
technique. In tapping mode, the microscope's cantilever 
tip oscillates at its resonant frequency as it is scanned 
across a surface. A feedback loop attempts to maintain 
the tip oscillation at a fixed amplitude by adjusting the 
average height of the tip above the sample surface; this 
average tip displacement is recorded as "height data." 
Transient changes in the oscillation amplitude, roughly 
a measure of the height derivative in the scan direction, 
are monitored as "amplitude data." The middle panels 
of Fig. 1 display height data, while the bottom panels 
display the corresponding amplitude data. 

The inevitable slow drifts of the cantilever tip during 
image acquisition, while leaving the amplitude data unaf- 
fected, can obscure shallow surface features in the height 
data. To compensate for these drifts, the height images 
in the middle panels of Fig. 1 were subjected to a high- 



2 





Sample A 




v= 9/2 




7/2 




1 1 



2.5 3.0 3.5 
B(T) 




► [110] 



FIG. 1. Top panels: Longitudinal resistance data taken 
at 50 mK in samples A and B. Dashed curves: resistance 
along [110]; solid curves: resistance along [110]. Strongly 
anisotropic transport is evident around v — 9/2. Middle pan- 
els: Height AFM images (16 x 16 (im) of the samples' parent 
wafers. The prominent mounds are 5 — 10 nm tall. Bot- 
tom panels: Corresponding amplitude AFM images. While 
for both samples the hard and easy transport directions at 
v — 9/2 lie along [110] and [110], respectively, the orienta- 
tions of their surface morphologies are perpendicular to one 
another. 



pass filter in the scan direction. Because this filtering 
risks a loss of information of surface height variations per- 
pendicular to the scan direction, four independent scans 
parallel to [lTO], [110], [010], and [100] were performed 
for each surface. These results (not shown here) confirm 
that the images of Fig. 1 faithfully represent all relevant 
morphology and that no spurious features related to the 
scan direction are present. Furthermore, the characteris- 
tic surface features revealed in Fig. 1 were found to cover 
each wafer's entire surface, indicating that their occur- 
rence is a generic property of the MBE growth process. 

Both surfaces in Fig. 1 exhibit shallow (< 10 nm), 
anisotropic surface features on the micron length scale. 
However, two important differences should be noted: 



first, the orientation of the surface features of the two 
samples are orthogonal to each other, and second, the 
general shapes of the surface features differ for the two 
samples. The surface of sample A, for example, is cov- 
ered by long and straight cigar-shaped mounds oriented 
with their major axes along [110]. Sample B's surface, on 
the other hand, exhibits more irregularly-shaped mounds 
clearly elongated along [110] and separated by narrow, 
meandering valleys. 

These results are typical of the variety of surface fea- 
tures seen on seven other wafers examined, where AFM 
images revealed varying degrees of mounding. Three 
of the seven showed mounds elongated along [110], two 
showed rounded mounds with virtually no preferred di- 
rection, and two showed mounds elongated along [110]. 
We emphasize, however, that in all the heterostructures 
the hard and easy axes of the resistance anisotropy in 
high LLs are along [110] and [110], respectively. Although 
the "quality" of the transport data does vary among the 
samples, it is sensibly correlated with the mobility of the 
2DES and shows no evident dependence on the surface 
morphology 

The presence of anisotropic mounding, or roughness, 
on thin GaAs epilayers grown by MBE on [001]-oriented 
GaAs substrates has been noted previously by Orme, et 
al. p4fl . The competition between step-flow growth and 
island nucleation, in conjunction with anisotropic surface 
diffusion, was found to result in micron-sized mounds 
elongated along [110]. Although Orme, et al. did not 
observe mounding oriented along [110], as we find in 4 of 
9 samples studied, we note that our heterostructures are 
much thicker and more complex (consisting of hundreds 
of layers of GaAs and AlGaAs) than simple GaAs epilay- 
ers. As Orme, et al. found the tendency to mound to be 
very sensitive to growth temperature, substrate miscut 
angle, and epilayer thickness, the greater diversity of our 
results seems less surprising. 

The lack of a consistent crystallographic orientation of 
the micron-scale surface roughness features on our sam- 
ples contrasts sharply with the highly consistent align- 
ment of the transport anisotropy axes exhibited by the 
2DES at v = 9/2, 11/2, etc., in the same samples. From 
this we conclude that micron-scale surface roughness does 
not reflect the symmetry breaking field responsible for 
the alignment of the anisotropic electronic phases in high 
LLs. Whether finer scale surface features or unseen sub- 
surface defects are involved remains to be determined. 
While we agree with Willett, et al. M|, that sufficiently 
severe surface roughness (whether as-grown or externally 
imposed) can influence 2DES transport, we do not find 
that the orientation of micron-scale surface features is al- 
ways correlated with transport anisotropics in high LLs. 

Another effect which might break the symmetry be- 
tween the [1T0] and [110] directions in a 2DES on a [001]- 
oriented surface is related to the shape of the potential 
well confining the electrons. Kroemer pi has noted that 
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while transport in bulk GaAs is invariant under it/ 2 ro- 
tations about the [001] axis, this may not be the case in 
a 2DES. If the potential confining the electrons to the 
[001] plane lacks inversion symmetry the transport coef- 
ficients need only be invariant under ir rotations. While 
it seems unlikely that the band structure effects result- 
ing from this lack of inversion symmetry could explain 
the large magnitude and sharp temperature and filling 
factor dependences of the high LL transport anisotropy, 
they may be sufficient to orient an otherwise sponta- 
neously generated anisotropic ground state. Rosenow, et 
al. [jlTj and Takhtamirov, et al. have stressed this 
point and have estimated that the reduced symmetry 
of typical single-interface heterojunctions leads to small 
(~ 0.1%) anisotropies in the 2D electron effective mass. 

The transport measurements of samples A and B in 
Fig. 1 allow for examination of this idea. Sample A con- 
sists of a single 30 nm GaAs quantum well embedded 
within Alo.24Gao.76As. Silicon doping layers positioned 
98 nm above and below the quantum well produce a sym- 
metrically confined 2DES. Sample B, on the other hand, 
contains a single GaAs/AlGaAs interface doped from one 
side only. The 2DES in this sample is therefore confined 
by an asymmetric, roughly triangular, potential. In spite 
of this structural difference the two samples show quali- 
tatively identical high LL transport features. In particu- 
lar, the orientation of the anisotropy axes relative to the 
crystal axes is the same in the two samples. While one 
might argue that uncontrolled microscopic irregularities 
would destroy the symmetry of square quantum wells, it 
has been convincingly demonstrated, in 2D hole systems, 
that the symmetry of the confinement potential can have 
readily measurable transport consequences under appro- 
priate circumstances |T^j2^]. Therefore, the lack of any 
significant differences in the high LL transport in sam- 
ples A and B suggests that the symmetry of the confine- 
ment potential is also not a major factor in determining 
the orientation of the underlying anisotropic electronic 
states. 

The transport data from sample A displayed in Fig. 
1, which were obtained without prior low temperature 
illumination, demonstrate that such illumination is not 
required to observe anisotropic transport in high Landau 
levels. Illumination is usually employed because it has 
been found empirically to often improve the quality of 
high field transport data. The precise way in which this 
occurs is not well understood, but in addition to typically 
increasing the 2D density and mobility, illumination also 
seems to improve the homogeneity of the electron gas. 
In any case, Fig. 2 demonstrates this improvement in 
sample A. The minima in the resistance measured along 
[110] at v = 9/2 and 11/2 are much deeper, the isotropic 
re-entrant integer quantized Hall states in the flanks of 
the Landau levels are much better developed |2p|,|2"l|] , and 
additional structure is evident. Nonetheless, these data 
show clearly that illumination is not an important factor 



in the origin or orientation of the anisotropic phases of 
2D electrons in high Landau levels. 
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FIG. 2. Effect of low-temperature illumination by a 
red light emitting diode on the transport characteristics of 
sample A at 50 mK. Top: Longitudinal resistance along 
[110] (dashed) and [110] (solid) before illumination in the re- 
gion 4 < v < 6. Bottom: The same measurement after illu- 
mination. Note the improvement in the degree of anisotropy 
at 1/2 filling and the increased prominence of the re-entrant 
integer states near 1/4 and 3/4 filling. 

A magnetic field component -By applied parallel to the 
plane of the 2DES provides an important tool for inves- 
tigating the magnitude of the native symmetry breaking 
mechanism. Such an in-plane field allows the controlled 
introduction of an additional symmetry breaker which is 
known to be capable of altering the high Landau level 
transport anisotropy axes Jll],|22|]. For example, only a 
relatively small Bu pointed along [110] is required to in- 
terchange the original principal axes of anisotropy, with 
the new easy and hard directions along [110] and [110] 
respectively. Prior to this interchange, there is a special 
value of the in-plane field, B^, where the transport in 
the sample becomes approximately isotropic. Remark- 
ably, if the in-plane field is instead directed along [110], 
the principal axes of anisotropic transport do not reverse. 

Hartree-Fock calculations within the undirectional 
CDW model can account for these observations: the cou- 
pling between Bu and the finite thickness of the 2D layer 
favors stripes aligned perpendicular to B\\ under typical 
circumstances |^3|,Q . Assuming that the high resistance 
direction is perpendicular to the stripes, this conclusion 
is in agreement with experiment. Similar calculations 
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have also been successful in explaining the more complex 
anisotropic transport behavior exhibited at high Bu by 
quasi-2D systems in wide quantum wells having two occu- 
pied subbands |^5| . These theoretical analyses also yield 
quantitative estimates of the field-induced anisotropy en- 
ergy, Ea, defined as the difference in energy (per elec- 
tron) of undirectional CDWs aligned parallel and per- 
pendicular to the in-plane magnetic field. By evaluating 
the anisotropy energy at Bp these calculations provide 
an estimate of the magnitude of the native symmetry 
breaking field which orients the transport anisotropy in 
the absense of an in-plane field. 

Figure 3 shows the effect of an in-plane magnetic 
field directed along [110] on transport in the vicinity of 
v = 9/2 in samples A and C (data from sample B has 
been published previously Jll|). These data, which are 
obtained by tilting the sample relative to the axis of a 
single superconducting solenoid, clearly reveal the inter- 
change of the anisotropy axes. For sample A a tilt of 
only 9 = 6° is sufficient to render the transport roughly 
isotropic. This corresponds to BT, ks 0.24 T at v = 9/2. 
By 9 = 10° the transport anisotropy is fully restored, al- 
beit with interchanged axes. Samples B and C exhibit 
switching fields of approximately BZ » 0.50 and 0.55 T 
respectively. 
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FIG. 3. Effect of a parallel magnetic field in the [110] di- 
rection on anisotropic transport at v — 9/2 of samples A (50 
mK) and C (25 mK). Tilting about the [110] axis results in 
a reversal of the principal axes of anisotropy. Dashed: longi- 
tudinal resistance along [110]. Solid: longitudinal resistance 
along [110]. 

Theoretical field-induced anisotropy energies are plot- 



ted in Figure 4. The many-body RPA/Hartree-Fock 
calculations are combined with a self-consistent local- 
density-approximation description of one-particle states 
p3j, obtained numerically for the specific structural pa- 
rameters of samples A, B, and C. In agreement with the 
experimental observations, the unidirectional CDWs pre- 
fer being aligned perpendicular to the in-plane magnetic 
field direction (Ea > 0). The agreement between ex- 
periment and theory on this point offers strong support 
to the overall CDW picture of the anisotropic phases in 
high Landau levels. After all, this theoretical result is 
not obvious a priori, as the interaction between a paral- 
lel magnetic field and the 2DES depends sensitively on 
the detailed finite- width profile of the 2DES and the na- 
ture of screening in high magnetic fields. Indeed, previ- 
ous calculations have shown that unlike strongly confined 
2D electron systems where the CDW orients perpendic- 
ular to Bh |23| ], in two-subband quantum wells the pre- 
ferred CDW orientation is parallel to the in-plane field 
|p5| . This observation helps explain the remarkable dif- 
ference in the dependence of Ea on Bu in samples A and 
B despite the samples having nearly identical 2D electron 
densities. In sample B, the unoccupied second subband 
is much closer to the chemical potential than in sam- 
ple A. Thus, being more similar to a system with two 
subbands occupied, a CDW in sample B should gain a 
smaller energy advantage by aligning perpendicular to a 
given in-plane field. 

The calculated field-induced anisotropy energies at the 
measured v — 9/2 switching fields BT< are 0.5, 0.2, and 
2.4 mK per electron in samples A, B, and C respectively. 
These energies are far smaller than the ~ 100 mK tem- 
perature scale that apparently governs the formation of 
the anisotropic electronic phase (i.e. the temperature at 
which transport becomes anisotropic). This disparity 
supports the model of a robust CDW oriented by a rel- 
atively weak native symmetry breaking field. Although 
the origin of the native symmetry breaking field remains 
unknown, these estimates of the anisotropy energy should 
offer useful constraints on proposed mechanisms. 

The experimental results reported here suggest that 
two possible symmetry breaking effects, anisotropic sur- 
face roughness and confinement potential geometry, are 
not important factors in determining the orientation of 
the newly discovered anisotropic phases of 2D electrons 
in high Landau levels. Nevertheless, other possibilities 
remain. For example, Fil |pq] has argued that the piezo- 
electricity of GaAs leads to lower energy for CDWs ap- 
proximately aligned along either [110] or [110]. Although 
this does not explain what further distinguishes between 
these two directions, it does establish a mechanism which 
breaks rotational symmetry and favors directions which 
do have relevance to experiment. 

In conclusion, we have addressed the still unresolved 
issue of how the anisotropic electronic states in high LLs 
become oriented relative to the crystallographic direc- 
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FIG. 4. Theoretical results of the anisotropy energy ver- 
sus an in-plane magnetic field for a unidirectional CDW at 
v = 9/2 in samples A, B, and C. The CDW periods in the 
three samples, at B\\ = 0, were calculated to be 108, 106, 
and 143 nm respectively. Evaluated at the switching field By , 
the anisotropy energy yields an estimate of the native pinning 
energy of a CDW phase. 



tions of the host semiconductor lattice. We have found 
that neither micron-scale surface roughness nor the sym- 
metry of the confinement potential are major factors. In 
conjunction with detailed theoretical calculations, tilted- 
field experiments support the overall charge density wave 
picture of the ground state of 2D electrons in high Lan- 
dau levels and yield quantitative estimates of the native 
anisotropy energy for three different samples. 
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